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ABSTRACT 

Two explanations exist for the short-lived radionuclides (Ti < 5 Ma) present 
in the solar system when the calcium-aluminum-rich inclusions (GAIs) hrst 
formed. They originated either from the ejecta of a supernova or by the in 
situ irradiation of nebular dust by energetic particles. With a half-life of only 
53 days, Beryllium-7 is then the key discriminant, since it can be made only by 
irradiation. Using the same irradiation model developed earlier by our group, we 
calculate the yield of ^Be. Within model uncertainties associated mainly with 
nuclear cross sections, we obtain agreement with the experimental value. More¬ 
over, if ^Be and ^“^Be have the same origin, the irradiation time must be short (a 
few to tens of years), and the proton flux must be of order F ~ 2 x 10^° cm“^ 
s“^. The x-wind model provides a natural astrophysical setting that gives the 
requisite conditions. In the same irradiation environment, ^®A1, ^®G1 and ®^Mn 
are also generated at the measured levels within model uncertainties, provided 
that irradiation occurs under conditions reminiscent of solar impulsive events 
(steep energy spectra and high ^He abundance). The decoupling of the ^®A1 and 
^°Be observed in some rare GAIs receives a quantitative explanation when rare 
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gradual events (shallow energy spectra and low ^He abundance) are considered. 
The yields of ^^Ca are compatible with an initial solar system value inferred from 
the measured initial "^^Ca/^°Ca ratio and an estimate of the thermal metamor¬ 
phism time (Young et ah 2005), alleviating the need for two-layer protoCAIs. 
Finally, we show that the presence of supernova-produced ®'^Fe in the solar ac¬ 
cretion disk does not necessarily mean that other short-lived radionuclides have 
a stellar origin. 

Subject headings: ^Be, irradiation, early solar system, cosmic rays, x-wind, short¬ 
lived radionuclides, meteorites 


1. Introduction 

An important result concerning the formation of the solar system is the discovery that 
radionuclides with half-lives of the order of a million years (Ma) and less were alive when 
Calcium-Aluminum-rich Inclusions (CAIs) and chondrules, the igneous components of chon¬ 
drites, were formed (see Russell, Gounelle & Hutchinson 2001 and McKeegan & Davis 2003 
for recent reviews). Furthermore, some short-lived radioactivities, ^°Be {Ti =1.5 Ma), ^®A1 
(0.74 Ma), ^®C1 (0.3 Ma), ^^Ca (0.1 Ma), ®°Fe (1.5 Ma) and possibly ®^Mn (3.7 Ma), ex¬ 
isted at levels (see Table 1) well above the expectations of models of continuous galactic 
nucleosynthesis (e.g. Meyer & Clayton 2000). These hndings have motivated two different 
kinds of explanations. The external model stipulates that the radionuclides were made in a 
supernova or in an asymptotic giant branch (AGB) star (e.g., Busso, Gallino & Wasserburg 
2003). In the case of a supernova, it has been proposed that the expelled material triggered 
the collapse of a nearby molecular cloud core, which led to the formation of the solar system 
(Cameron & Truran 1977). More recently, an Orion-like environment has been suggested 
for the birth of the solar system (Hester et al. 2004). The internal irradiation model con¬ 
jectures that the radionuclides were created by the irradiation of the solar environment by 
energetic protons, ^He and ^He nuclei, accelerated by gradual or impulsive events near an 
active young Sun (Lee et al. 1998; Goswami et al. 2001; Gounelle et al. 2001; Marhas et 
al. 2002; Leya et al. 2003). Thus, establishing the origin of the short-lived radionuclides 
can provide important constraints on the astrophysical setting of the Sun’s birth, stellar 
nucleosynthesis models, irradiation processes around young stellar objects, and early solar 
system chronology (Gounelle et al. 2005a,b; Kita et al. 2005). 

Berylhum-7 decays to ^Li with a half-life of 53 days. Because of this very short half- 
life, it must have an internal irradiation origin, since any ^Be produced outside the nascent 
solar system would have decayed long before being incorporated into the hrst condensing 
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solids. While the mere discovery of ^Be would be dehnite proof for in-situ irradiation, 
measurement of its abundance, together with that of other irradiation products such as ^'^Be 
(McKeegan Chaussidon and Robert 2000), would provide important information concerning 
the irradiation conditions (duration as well as flux, energy spectra and isotopic composition 
of the accelerated particles). 

The experimental detection of ^Be has proven to be a long and difficult quest, largely be¬ 
cause of the high mobility of its daughter element, Li. A hint of the presence of ^Be was first 
seen in an Allende CAI (USNM 3515) at a high relative abundance, ^Be/®Be ~ 0.1 (Chaus¬ 
sidon, Robert, & McKeegan 2002). However, a dehnitive isochron could not be established 
for this CAI, which has been highly disturbed by secondary events. Recently, Chaussidon, 
Robert & McKeegan (2004, 2005) identihed 37 spots within Allende CAI 3529-41 for which 
the Li isotopic record has been undisturbed since crystallization. They determined a well- 
dehned isochron for these 37 points and established that ^Be was alive in the early solar 
system at a level ^Be/®Be = (6.1 ± 1.3) x 10“^ {2a error). 

The goal of this paper is to examine the yields of ^Be in the context of the irradiation 
model we have developed over the past several years (Shu et al. 1996; Shu et al. 1997; Lee 
et al. 1998; Gounelle et al. 2001). We will establish additional constraints on irradiation 
models imposed by the presence of ^Be and ^'^Be in the early solar system. Calculations for 
the production of ^®C1, a newly discovered radionuclide (Lin et al. 2005) with a half-life of 
0.3 Ma, are also presented, as well as new data and updated computations for ^®A1, ^^Ca 
and ^^Mn. We also discuss the implications of the high abundance of ®°Fe in the early solar 
system. 


2. The Irradiation Model 
2.1. Outline of the Model 

We use the irradiation model in the form developed by Gounelle et al. (2001) to calculate 
the yield of ^Be. Their preferred case (c.f. Figure 2d of that paper) is adopted, and the 
parameters are kept the same except for the spectral index p and the ^He/^H ratio, which 
are varied within reasonable ranges. Case 2d refers to the geometry of the assumed core¬ 
mantle structure of protoCAIs: the minimum and maximum core sizes vary between 0.0050 
and 2.5 cm, as natural CAIs do, and the mantle thickness is fixed at 0.28 cm (Shu et al. 2001). 
Here we provide the main characteristics of our model and discuss how it is implemented for 
the new calculations of ^Be and ^®C1. The very short half-life of ^Be (53 d vs ~ Ma for other 
radionuclides) calls for a special treatment. Additional details on the irradiation model may 
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be found in Gounelle et al. (2001). 

Of basic importance in irradiation models are the absolute yields of the short-lived 
radionuclides. By absolute yield we mean the abundance of a radionuclide calculated relative 
to a hxed, stable-isotope reference (e.g. ^^Al/^'^Al), not abundances of radionuclides relative 
to one another. The absolute yields depend linearly on the product of the number flux of 
cosmic rays, TcR) and the irradiation time, At. These quantities can vary over many orders 
of magnitude, depending on the astrophysical environment, implying changes in the short- 
lived-radionuclide yields of many orders of magnitude. Other parameters, like the composite 
structure of the targets, the shielding, or the cosmic ray properties, induce variations in the 
yields of the short-lived radionuclides that are no larger than a few orders of magnitude (Lee 
et al. 1998; Goswami et al. 2001; Gounelle et al. 2001; Marhas et al. 2002; Leya et al. 
2003). 

In our model, the number flux of protons is scaled to X-ray observations of young stellar 
objects (Lee et al. 1998), with Fp {E > 10 MeV) ~ 1.9 x 10^° cm^ s“L The irradiation of 
proto-GAIs occurs within the reconnection ring (Lee et al. 1998), and the effective number of 
proton passages in the reconnection ring is Af = 1.6. The irradiation time depends linearly 
on the size of the proto-GAIs. With a typical irradiation time of 20 yr for a 1 cm protoGAI, 
the irradiation time varies between a few years and a few tens of years, since protoGAIs vary 
between 0.2805 and 2.78 cm (Gounelle et al. 2001; Shu et al. 2001). Homogenization of 
the irradiation products is assumed to take place through frequent evaporation-condensation 
episodes that, in a steady-state picture, erase the calculated differences of the yields between 
small and large GAIs. The hnal yield of a short-lived radionuclide is obtained by integrating 
the irradiation products over the size distribution of the total population of GAIs, reflecting 
the continuous homogenization of irradiation products made possible by the high frequency 
of large flares compared to the irradiation duration (Grosso et al. 1998; Wolk et al. 2005). 

The differential distribution of protons, N{E), is assumed to be a power-law in energy. 
The spectral index p in N{E) oc E~'p has a range between 2.7 and 5. The number of 
^He nuclei relative to protons is hxed at 0.1, but the ^He/^H ratio varies between 0 and 1. 
Impulsive events are characterized by relatively steep proton spectra (higher p) and large 
abundances of ^He, while gradual events have shallower proton spectra (lower p) and smaller 
abundances of ^He (Reames 1995). We envision two different structures for proto-GAIs, a 
core-mantle conhguration and a homogeneous chondritic composition (Gounelle et al. 2001). 

For the production of ^®G1, we consider proton, alpha and ^He induced reactions on a 
diversity of targets such as Gl, S and K. The adopted average chondritic abundance of G1 is 
704 ppm (Lodders 2003) while the adopted median GAI abundance is 390 ppm (Sylvester et 
al. 1993). 
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For the production of ^Be (and ^'^Be) we consider proton, alpha and ^He induced re¬ 
actions on Because its mean-life r (76 days) is short compared to typical irradiation 
times At, ranging from a few years to a few tens of years, ^Be yields need to be calculated 
by balancing spontaneous decay against production by spallation reactions. Thus, for ^Be, 
and any short-lived radionuclide for which r At, the yield is the production rate times r 
rather than times At, which applies when r ^ At. 


2.2. Updated Cross Sections 

The main new input data are the production cross sections shown in Figure 1 and 
Figure 2. The ^Be cross sections are based on numerical simulations using fragmentation and 
Hauser-Feshbach codes and on experimental data points (indicated by circles) from Landolt 
& Bornstein (1994). The ^‘^Be-production cross-sections have been updated since the work 
of Gounelle et ah (2001). The ^®0(p,x)^°Be cross section is from Sisterson et ah (1997); the 
^®0("^He,x)^°Be cross section is from the experimental study of Lange et al. (1995) that had 
been previously overlooked; and the ^®0(^He,x)^°Be cross section is from a new simulation. 

The ^^Cl-producing cross sections (Figure 2) are based on numerical simulations using 
fragmentation and Hauser-Feshbach codes. In the absence of experimental data, the input 
model parameters have been chosen so that the calculations reproduce data obtained for lower 
and higher masses (e.g. Al. Ca, Fe) where some data exist. Because of the dependence of the 
cross sections on the structure of the nuclei involved, it is difficult to give precise uncertainties 
for these cross sections. From comparisons with experimental data for other reactions, we 
estimate an uncertainty of 50% for the maximum of the cross section and a factor of 2 for 
cross sections that are changing significantly with energy. 

Most of the ^He cross-sections used by Lee et al. (1998) and by Gounelle et al. (2001) 
for the production of short-lived radionuclides came from theoretical nuclear physics codes 
and not from direct laboratory measurements. To help reduce the uncertainties based on 
this over-reliance on theory, Fitoussi et al. (2004) measured the cross section for the reaction 
^"^Mg(^He,p)^®Al that is the main source of ^®A1. In addition to using this experimentally 
determined cross section, we also include proton and alpha reactions for the production of 
^®A1 and ^^Ga that were previously ignored (for impulsive events) on the basis that their 
contributions to the ^®A1 and ^^Ga yields were small compared with ^He. These cross sec¬ 
tions are: ^®Mg(p,n)^®Al, ^'^Al(p,pn)^®Al, ^®Si(p,2pn)^®Al, ^^Mg(Q;,pn)^®Al, "^^Ga(p,pn)^^Ga, 
^°Ga(Q;,^He)"^^Ga, ^°Ga(Q;,^H)^^Sc, and they can be found in Lee et al. (1998). 

To summarize, we now use an improved and a more complete set of cross sections 
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compared to Gounelle et al. (2001). Most of the key cross sections for the production of ^Be, 
^°Be and are now experimentally constrained. Among the cross sections important for 
our calculations, only ^®0(^He,x)^°Be, ^°Ca(^He,pn)^^Ca as well as the ^®C1 producing cross 
sections lack experimental data. The uncertainties in the cross sections in Figure 1 (where 
experimental points are labeled with hlled circles) are a factor of the order of a few, and 
possibly much larger for cross sections without experimental measurements such as the ones 
presented in Figure 2. 


2.3. Results of the Calculations 

Table 2 gives the calculated absolute yields of ^Be as a function of the spectral index p 
and the ^He/^H ratio for two different cases: a core-mantle and a homogeneous chondritic 
compositional structure. The yields have been normalized to the experimental value obtained 
by Chaussidon, Robert and McKeegan (2005). Figure 3 shows the (relative) ^°Be/^Be ratio 
as a function of the ^He/^H ratio for a range of values of the spectral index p. Using the 
absolute ^Be yields in Table 2, this is equivalent to giving the absolute yield of ^'^Be. The 
result obtained with the same parameters as Gounelle et al. (2001), corresponding to ^He/^H 
= 0.3 and p = 4, is marked by a hlled circle in Figure 3. 

Figure 4 shows the irradiation yields for the short-lived radionuclides of primary interest, 
^Be, ^'^Be, ^®A1, ^®Gl,"^^Ga, ^^Mn, for the preferred parameters, p = 4 and ^He/^H = 0.3, and 
for the two different target compositions (core-mantle and chondritic). Figure 5 presents the 
results for a variety of spectral parameters. The case p = 2.7 and ^He/^H = 0 represents 
typical gradual hares, while the case p = 4 and p = 5 with ^He/^H = 0.3 represents impulsive 
hares. We now discuss the implications of these results. 


3. Discussion 

3.1. The Irradiation Origin of ^Be 

The detection of ^Be, with a half-life of 53 d, in meteorites (Ghaussidon, Robert and 
McKeegan 2005) provides dehnitive evidence for in situ irradiation of pre-meteoritic solids 
and for local irradiation models. In fact, an earlier theoretical calculation, using preliminary 
cross section data (Gounelle et al. 2003), predieted the experimental result (Ghaussidon, 
Robert & McKeegan 2004) to within a factor of two. An explanation of ^Be in terms of 
some other process (e.g., Gameron 2003) is possible only if GAIs themselves formed outside 
the solar system and were never subsequently re-heated to reset the beryllium-lithium clock. 
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In such a case, however, the entire dating of the solar system, based on the nearly identical 
radiochemical ages of all CAIs, would be called into question. 

For reasonable cosmic ray parameters {p ranging between 2.7 and 5, and ^He/^H ranging 
from 0 to 1), the ^Be yields in Table 2 vary between 0.12 and 5.6 times the experimental 
value. The preferred model of Gounelle et ah (2001), with p = 4 and ^He/^H = 0.3, gives 
^Be/®Be = 7.5 x 10“^, which compares favorably with the measured value of ^Be/®Be = 6.1 
X 10“^. As can be seen in Table 2, the ^Be yields do not depend very much on the chemical 
composition of the target (usually no more than a few per cent). This is because they depend 
on the 0/Be ratio, and refractory Be undergoes minimal chemical fractionation. Although 
the agreement between the model and the experimental value is very good, uncertainties in 
the model, especially in the nuclear cross sections (see section 2), call for caution. 

A possible source for variations in the yields of ^Be are fluctuations of the energetic 
particle flux. Because the relevant timescale for ^Be synthesis is tens of days, ^Be is more 
likely to be sensitive to such fluctuations than, say, ^°Be with a half-life of 1.5Ma. Signihcant 
X-ray flares are observed for Sun-like young stellar objects with ages in the 1-2 Ma range on 
the order of several days (Wolk et al. 2005). Thus '^Be synthesis in natural samples is likely 
to be more variable than other short-lived radionuclides. 

A basic reason for the good agreement between the present model calculations and 
experiment is that the production of ^Be requires a high accelerated particle flux. As pointed 
out in §2), the yield of ^Be with its short half-life does not depend on the irradiation time 
At but on the mean life r. It also depends only moderately on the cosmic-ray parameters, 
as can be seen in Table 2. Consistent with the results in Table 2, the only way to produce 
^Be at a signihcant level is to have a high hux of accelerated particles. The value adopted by 
Gounelle et al. (2001) and used here, Fp{E > Eiq) ~ 1.9 x 10^° cm^ s“^, is consistent with 
the measured X-ray properties of Sun-like stars (e.g., Wolk et al 2005). Our local irradiation 
scenario is also naturally conducive to high accelerated-particle huxes because the particles 
are conhned to as well as produced within the reconnection ring (the irradiation zone) by 
strong magnetic helds (see 3.5). 


3.2. The Case of ^°Be 

Our results for ^'^Be are shown in Figures 3, 4 and 5. The calculated ^°Be/^Be ratio 
depends moderately strongly on the cosmic-ray parameters (index p and ^He/^H), when 
^He/^H > 0.01, as shown in Figure 3. The preferred model of Gounelle et al. (2001) gives a 
satisfactory account of the experimental ratio: (^°Be/'^Be)model = 0.88 vs. (^°Be/^Be)exp = 



0.16. The over-production of ^°Be by a factor of ~ 5 is not that large, given the uncertainties 
in the cross sections. Moreover, it is reduced if we use a somewhat steeper spectral index, 
e.g., p = 5 as shown in Figure 5. Alternatively, the ^°Be yield would be lower if the irradiation 
time was lower. However, this modihcation would also decrease the yields of other short-lived 
radionuclides such as ^®A1. 

The ratio ^'^Be/^Be does not depend at all on the number flux (see section §2.1) and 
very little on the target composition (see Table 2). It depends mostly on the irradiation time, 
which is estimated in our 2001 model to be a few years to a few tens of years. The irradiation 
time is essentially the residence time for proto-CAIs of a certain size in the reconnection ring 
before ejection to asteroidal distances by the fluctuating X-wind (Shu et ah 2001). Such time 
scales are also suggested by observations of the intervals between knot-forming episodes in 
the optical jets of young stellar objects (Reipurth and Bally 2001). We cannot emphasize too 
strongly that the resulting mechanical model (Shu et ah 2001) yields in order-of-magnitude 
exactly the irradiation time to reproduce the ^°Be/^Be ratio. 

Desch, Connolly and Srinivasan (2004) have challenged the idea that ^‘^Be was formed in 
the solar system by irradiation. They argue instead that most of the solar system ^°Be orig¬ 
inated as ^'^Be nuclei in galactic cosmic rays stopped in the progenitor molecular cloud core. 
This conclusion depends sensitively on a number of assumptions made in their calculations, 
e.g., in the variation of the mass-column density S(ci7,t) in a collapsing molecular-cloud 
core, taken from an earlier numerical study by Desch & Mouschovias (2001). Desch et al. 
(2004) regard S as being spatially uniform, whereas Desch & Mouschovias (2001) found S 
to decrease considerably with increasing distance w from the center of the cloud core. The 
core has a mass (45 M©), signihcantly in excess of what is expected for the formation of a 
solar-mass star. Furthermore, the time scale for core formation and ^°Be accumulation is 
~ 10 Ma, an order of magnitude longer than allowed by the statistics of nearby molecular 
cloud cores with and without embedded stars (e.g., Lee & Myers 1999). Even so, Desch et 
al. (2004) still need to increase the galactic cosmic ray flux relevant to present values by an 
ad hoc factor of 2 to produce enough ^°Be. The long evolutionary time also leaves no room 
for the injection into the nascent protosolar cloud of other short-lived radionuclides like ^®A1, 
^^Ca and ^^Mn from a stellar event. Because of its ad hoc nature, the proposal of Desch, 
Connolly and Srinivasan (2004) has few advantages over the more common view that ^'^Be 
has an in-situ irradiation origin (e.g., McKeegan, Chaussidon & Robert 2000; Gonnelle et 
al. 2001; Marhas et al. 2002; Huss 2004). 



9 


3.3. The Co-production of ^®A1, ^®C1, and ®^Mn 

Our preferred model (case 2d of Gounelle et al. 2001) with p = A and ^He/^H = 0.3) 
can reproduce within a factor of a few the measured solar system abundance of ^Be, ^°Be, 
26Ai^ 36 qi^ ^^Ca and ^^Mn (the thick solid curve in Figure 5). Chlorine-36 is slightly under 
produced relative to its experimental value but, given the uncertainties of the model, the 
agreement is satisfactory. The adoption of a new experimental cross section for the reaction 
^"^Mg(^He,p)^®Al does not modify the yield of ^®A1 by more than a factor of 2, keeping the 
production of this radio-isotope in line with what had been calculated using theoretical cross 
sections. Taking into account the previously ignored proton and ^He induced cross sections 
for impulsive flares does not change the yields by more than ten percent. However, the 
adoption of an experimental cross section for ^®0(^He,x)^‘^Be has slightly increased the ^'^Be 
yield, degrading the quality of the £t with the experimental data. Adoption of a somewhat 
steeper spectral index (p = 5 instead of p = 4, all other parameters remaining the same) 
would put ^°Be and ^®C1 in line with the measured solar system value. This would not change 
much the yields of the other short-lived radionuclides (Figure 5), except for ^^Ca discussed 
below. 

The yields of all of the main short-lived radionuclides (including the newly calculated 
chlorine-36) are independent of composition except for ^^Ca, as may be seen in Figure 4. 
The low "^^Ca measured value, "^^Ca/^°Ca = 1.5 x 10“® (Srinivasan et al. 1994), can be 
reproduced by our model only when a core-mantle chemistry is adopted (Shu et al. 1997; 
Gounelle et al. 2001). However, there are hints that the solar system initial "^^Ga/^°Ga may 
have been considerably higher (McKeegan et al. 2004) than the presently measured value. 

Recent experimental findings suggest that the initial ^®A1/^'^A1 ratio of the solar system 
was higher than the measnred ’’canonical” valne 4.5 x 10“^ (MacPherson et al. 1995), with 
an initial ^®A1/^'^A1 of the order ~ 7 x 10“^ (Yonng et al. 2002; Galy et al. 2004; Yonng 
et al. 2004). Snbseqnent redistribntion of Mg isotopes within the GAI would then lead to 
an apparent initial ratio of 4.5 x 10“®, i.e., the canonical valne. Given the nncertainties of 
onr irradiation model, an increase of ~ 20 % in the initial ^®A1/^^A1 ratio would not change 
the overall level of agreement between the calculations and the measnrements. However, it 
might partly solve the overprodnction of ^^Ga observed when a pure chondritic composition 
is adopted (Figure 5). 

Galcium-41 decays into (see Table 2). If the K isotopes had been subjected to 
the same isotopic resetting as the Mg isotopes, this would mean that the initial "^^Ga/^°Ga 

t(26A0 t(26A0 

ratio might have been higher by a factor of (^^Alinitiai/^^Aheset) = (7/4.5) = 26, 

where we have simply applied exponential decay law to both ^®A1 and ^^Ga. It is therefore 
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very likely that the initial ^^Ca/‘^°Ca ratio was closer to 3.9 x 10“’^ than to 1.5 x 10“®. 
Such an early solar-system value would be consistent, within errors, with our computed 
irradiation yield of ^^Ca without the ad hoc adoption of a core-mantle structure (see Figure 
6). Isotopic redistribution would be less important for other short-lived radionuclides because 
they have longer half-lives than ^®A1. Chlorine-36 could also have been higher by a factor 
of ~ 3.8 because of this thermal metamorphism, slightly degrading the good fit of Figures 
5 and 6. Chaussidon, Robert and McKeegan (2005) have observed such a post-magmatic 
redistribution for Li isotopes, but have been able to deconvolve its effect from the ^Be decay. 

3.4. Decoupling of A1 and 

Mahras, Goswami & Davis (2002) have detected the past presence of ^°Be in some rare 
refractory inclusions containing calcium aluminum oxides known as hibonites in which ^®A1 
and "^^Ca (Sahjipal et ah 1998) are conspicuously absent. They interpret this “decoupling” 
between ^'^Be and ^®A1 (and ^^Ca) as contradicting irradiation models that simultaneously 
produce ^'^Be, ^®A1 and "^^Ca. 

Hibonites from CM2 chondrites, which show this decoupling between ^'^Be and ^®A1, 
share with FUN inclusions the property of bearing large isotopic anomalies in neutron-rich 
isotopes. Such FUN inclusions also formed with little initial ^®A1, a fact accounted for by 
Lee et al. (1998) who suggested that FUN CAIs never experienced the impulsive-flare events 
that are rich in ^He before being ejected by the X-wind. However, CAIs that never enter the 
reconnection ring can still be exposed to gradual-flare events that occur at larger distances. 
Such gradual events are capable of producing large amounts of ^°Be without co-producing 
^®A1 (or "^^Ca), according to the results shown in Figure 5. Thus, the occurrence of rare 
CAIs containing abundant ^'^Be but little ^®A1 or "^^Ca may have an explanation in exposure 
to gradual but not impulsive events.We note that the anomalies in stable isotopes are not 
produced by irradiation and are probably of nucleosynthetic origin (Fahey et al. 1987a). 

The recent measurement of ^‘^Be initial value in the hibonite HAL makes our hypothesis 
amenable to confrontation with experimental data. Marhas and Goswami (2003) have mea¬ 
sured ^°Be/®Be of (8.0 ± 3.3) x 10“"^ in HAL. Using ^^Al/^'^Al = (5.4 ± 1.3) x 10“® (Fahey 
et al. 1987b), we get = (8.1 ± 4) x 10^ for HAL. For gradual flares (p = 2.7 and 

^He/^H = 0), we obtain = 4.7 x 10^. The agreement within errors between the 

measured and the calculated ratio of ratios, indicates that our model can success¬ 

fully account for the decoupling between ^®A1 and ^'^Be measured by Marhas, Goswami and 
Davis (2002). Our model predicts that the ^Be/^'^Be ratio should be lower than normal in 
hibonites and FUN inclusions. For the sake of simplicity in the actual calculations, we have 
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assumed the same scaling of particle fluxes from X-ray observations for both gradual and 
impulsive events, but they could well be different. In other words, one should not deduce 
from Figure 5 that the absolute ^'^Be abundance is expected to be higher in gradual events 
compared to impulsive events, or that the ^'^Be/^Be should be higher in FUN inclusions than 
in normal inclusions. 


3.5. The Astrophysical Environment 

Irradiation in the reconnection ring produces the relatively high energetic particles fluxes 
that characterize our model. It is in this region that the magnetic helds generated by the 
young star and the accretion disk reconnect and accelerate protons and helium nuclei to 
energies of order of lOMeV or more. The reconnection ring is conhned by magnetic fields 
that prevent the ionized particles from escaping. The short irradiation times used in the 
model correspond to the residence times of proto-CAIs calculated for the reconnection ring 
(Shu et al. 2001), and are in harmony with the time scales observed for jets from active young 
stellar objects. The absolute proton flux is obtained from the observed X-ray luminosities of 
young stellar objects using a normalization based on measurements of X-rays and particle 
fluxes for the contemporary active Sun. The X-ray luminosity of active young stellar objects 
estimated by Lee et al. (1998), and used by Gounelle et al. (2001) as well as in this paper, is 
supported by recent measurements of the X-rays emitted by young stellar objects in the Orion 
Nebula Cluster (ONC) by the Chandra satellite (Feigelson, Garmire, & Pravdo 2002; Wolk 
et al. 2005). The proton flux used in our irradiation calculations for the reconnection ring 
(cited above as Fp{E > lOMeV) ~ 1.9 x 10^° cm^ s“^), corresponds to an X-ray luminosity 
of ~ 5 X 10^° erg s“^ (Lee et al. 1998; Gounelle et al. 2001). This is in line with the median 
value for flares observed by Wolk et al. (2005), Lx ~ 10^^ erg s“^ in an unprecedented Ms 
exposure by Chandra oi the ONC. These flares occur on an average repetition time of several 
days. Much more luminous flares are observed on longer time scales of the order of one year. 
Of course these observations of flares in the ONC refer mainly to T Tauri stars with ages 
~ 1 — 2 Ma, whereas much of the irradiation most likely occurs earlier during a much more 
active stage of accretion. Following the discussion of Lee et al (1998), the flares in that epoch 
were probably even more powerful. 

The irradiation model developed by Goswami, Marhas and Sahjipal (1998), and ex¬ 
tended by Marhas, Goswami and Davis (2002), does not feature a favorable environment 
for the high particle fluxes needed to explain the Be isotopes. These authors choose the 
asteroid belt, 2-3 AU distant from the young Sun, as the place to irradiate meteoritic precur¬ 
sors by protons and alpha particles. Because they consider the rare hibonites that contain 
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irradiation-produced and no to be the first solar system solids to form (Sahjipal 
and Goswami 1998), irradiation in their model takes place well before the accretion disk 
was dissipated. Because a lOMeV proton is stopped when it has encountered 0.2 grcm“^ of 
matter (Reedy 1990), the flux at asteroidal distances is expected to be well below the level 
needed to form ^Be and other short-lived radionuclides, especially when the inverse-square 
dilution of the stellar energetic particles is taken into account. The irradiation model devel¬ 
oped by Leya, Halliday and Wieler (2003) adopts the x-wind model as a general framework 
for irradiation close to the young Sun, followed by wind-borne delivery of the irradiated ma¬ 
terials to the asteroid belt. However, Leya et al. (2003) do not adopt a high energy particle 
flux or an irradiation time from observations or first principles, but calculate yields relative 
to the ^‘^Be/^Be ratio. 


3.6. The Stellar Origin of ®°Fe 

Neutron-rich ®°Fe is almost impossible to synthesize in the internal irradiation scenario 
because of the low abundance of appropriate (neutron-rich) projectiles and targets (Lee et 
al. 1998). However, it is well known that Type-H supernovae can produce large amounts 
of ®°Fe (Timmes et al. 1995; Rauscher et al. 2002). Using updated nuclear reaction yields 
and improved stellar physics, Rauscher et al (2002) estimate that Type-H supernovae can 
synthesize between 2.4 x 10“^ and 16 x 10“^ Mq of ®'^Fe, depending on the mass of the 
progenitor, somewhat higher than Timmes et al. (1995), 0.29-10.5 x 10“^ M©. 

The initial solar system value of ®°Fe/^®Fe is relatively well constrained since it has now 
been identihed in the silicate portion of chondrules that were presumably formed early in 
the history of the solar system (Tachibana & Huss 2003). Tachibana et al. (2005) estimate 
that the initial ®'^Fe/^®Fe was 5-10 x 10“’^. Smaller previous estimates were inferred by 
Shukolyukov & Lugmair (1996), ®°Fe/®®Fe = 6 x 10“®, from ^‘^Fe identihed in planetary 
differentiates that formed relatively late. The Tachibana et al. (2005) value is compatible 
with the initial ^‘^Fe/^^Fe = 9.2 x 10“® estimated by Mostefaoui et al. (2005) from measure¬ 
ments in Semarkona troilite. Using a solar system abundance of ^®Fe, ^®Fe/^H = 3.16 x 10“® 
(Fodders 2003), the initial ®°Fe content of the solar system lies between 0.95 x 10“® M© and 
1.9 X 10-9 Mq. 

The updated ^^Fe yields for Type-H supernovae and the new initial ^^Fe content of 
the solar system have important implications for the origin of ^^Fe and other short-lived 
radionuclides. They indicate that a supernova can deliver enough ^^Fe to the solar system’s 
progenitor molecular cloud core to account for the measured value. Other short-lived ra¬ 
dionuclides are not equally easily co-delivered at their initial abundances. The contribution 



13 


of a supernova, or any star, to the inventory of short-lived radionuclides in the early solar 
system can be characterized by a number of factors: the stellar yields, the mass fraction / 
of ejecta that is injected in the protosolar accretion disk, and the time interval A between 
nucleosynthesis and delivery. For example, consider the model SN25P of Rauscher et ah 
(2002) that produces 1.6 xl0“"^ Mq of ^‘^Fe by adopting an injection fraction / = 1.4 x 
10“"^ (2.8 X 10“^) and a time interval A = 7 Ma. By construction, such a supernova injects 
0.95 X 10“® Mq (1.9 X 10“^ Mq) of ®°Fe in the protosolar accretion disk, in line with the 
observed value. Such a supernova would also deliver ^®A1, ^®C1 and "^^Ca at levels well below 
the estimated abundances for the early solar system. 

No stellar model yet published can inject into a model protosolar accretion disk the whole 
inventory of short-lived radionuclides (^®A1, ^®C1, ^^Ca, ®°Fe). In his study of a specihc 25 Mq 
supernova model that delivers ^®A1, ^^Ca and ^'^Fe at the correct abundance, Meyer (2005) 
fails to deliver ^®C1 by two orders of magnitude. Moreover, even the highest-mass stars 
require at least two million years of normal stellar evolution before they supernova, and low- 
mass molecular cloud-cores that form sunlike stars in Orion-like environments (e.g., Hester 
et ah 2004) have disappeared before this time. Similarly, although specihc nucleosynthetic 
models of AGB stars (M = 1.5 Mq and Z = Zq/Q) have the potential to co-deliver ^®A1, 
^^Ca and ®°Fe to the solar system, they fail to give ^®C1 by two orders of magnitude and do 
not produce any ^^Mn (Gallino et ah 2004). Moreover, the association of a molecular cloud 
core and an AGB star is an improbable astrophysical event (e.g., Kastner and Myers 1994). 
Gontrary to what has been often argued (Zinner 2003; Hester et ah 2004), the presence of 
®°Fe at a high level in the nascent solar system does not necessarily imply that ^®A1 and other 
short-lived radionuclides also have a stellar origin. A wide range of supernovae, distant in 
space and time from the early solar system, can deliver high levels of ^‘^Fe without producing 
signihcant amounts of other short-lived radionuclides. 


4. Conclusions 

The discovery of ^Be (half-life 53 days) has important implications for early solar-system 
processes. First it dehnitely establishes that irradiation took place at an early stage in its 
formation. Second, it requires high accelerated particle fluxes, Fp (E > 10 MeV) ~ 2 x 
10^° cm^ s“^. When the measurements of ^Be are combined with those for ^“^Be, short 
irradiation times of the order of a few years to a few tens of years are adduced. The 
model developed by Lee et ah (1998) and Gounelle et al. (2001) satisfy these demanding 
conditions for high accelerated particles fluxes and short irradiation times. Our irradiation 
model produces reasonable abundances for other short-lived radionuclides, such as ^®A1, 
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36C1, and ^^Mn, as well as the Be isotopes. The overall scenario is also consistent 

with the evidence for energetic photo-irradiation of early solar system rocks from the mass- 
independent fractionation observed in stable sulfur and oxygen isotopes (Rai, Jackson, & 
Thiemens 2005). 

The decoupling between ^‘^Be and ^®A1 observed in a few rare CAIs can be accounted for 
by rare gradual events. For "^^Ca, a core-mantle chemistry is required if the initial value is as 
low as ^^Ca/^°Ca = 1.5 x 10“®. If the initial abundance of "^^Ca were a factor of 25 higher, 
as suggested by McKeegan et ah (2004), its abundance could be accounted for, within the 
model uncertainties, by a homogeneous chondritic composition, rather than a core-mantle 
structure. As noted earlier by Lee et al. (1998), local irradiation models fail to produce ®°Fe 
at the measured early solar system level, ^‘^Fe/^^Fe = 5-10 x 10“^ (Tachibana et al. 2005). 
A distant supernova is a likely origin for ®°Fe, as previously suggested by Lee et al. (1998). 
We note that recent updated models of supernova nucleosynthesis produce large amounts of 
®°Fe relative to other short-lived radionuclides, implying that it is possible for supernovae to 
inject ®°Fe in the solar system without co-delivering other short-lived radionuclides. 
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Fig. 1.— Adopted nuclear cross sections for the reactions ^®0(CR,x)^Be and ^®0(CR,x)^'^Be. 
The solid curves are the ^Be cross sections, and the dashed lines are the ^°Be cross sections. 
The incident projectiles, CR = p, ^He, and a are labelled on the graph. The nuclear cross 
sections are a combination of experimental data and Hauser-Feshbach fragmentation codes. 
Experimental data points are indicated by full black circles. 
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Fig. 2.— Adopted nuclear cross sections for the proton (thick line), a (dotted line) and ^He 
(dashed line) induced reactions for the production of ^®C1. The cross sections are based on 
numerical simulations using fragmentation and Hauser-Feshbach codes. In the absence of 
experimental data points, the input parameters have been chosen so that the calculations 
reproduce data obtained for lower and higher masses (e.g. Al. Ca, Fe) where some data 
exist. 





Fig. 3.— Plot of the ^‘^Be/^Be ratio, normalized to the measured solar system value, as a 
function of the spectral parameter p and the ^He/^H abundance ratio. The horizontal dotted 
line denotes the experimental value calculated from Chaussidon, Robert and McKeegan 
(2005) and McKeegan, Chaussidon and Robert (2001). The hlled circle corresponds to the 
preferred model of Gounelle et al. (2001), i.e with p = 4 and ^He/^H = 0.3. As in the case 
of ^Be (Table 2), the dependence on the target composition is weak and is not shown here. 
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Fig. 4.— Yields of ^Be, ^°Be, ^®A1, ^®C1, and for fixed spectral parameters (p 

= 4, ^He/^H = 0.3) and two different chemical compositions (core-mantle and chondritic). 
The yields are normalized to the experimental valnes (see Table 1), corresponding to the 
horizontal dotted line. 
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7Be 10Be 26AI 36CI 41 Ca 53Mn 

Fig. 5.— Yields of ^Be, ^®A1, ^®C1, and for fixed core-mantle structure (case 

2d of Gounelle et al. (2001)) and a variety of spectral parameters. Impulsive flares have 
p = 4 (thick line) or p = 5 (thin line), and ^He/^H = 0.3 while gradual flares have p = 2.7 
and ^He/^H = 0 (dashed line). The yields are normalized to the experimental values (see 
Table 1), corresponding to the horizontal dotted line. 
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Fig. 6.— Yields of ^Be, ^°Be, ^®A1, ^®C1, and for a chondritic chemistry and the 
spectral parameters p = 4 and ^He/^H = 0.3. The yields are normalized to the experimental 
valnes (see Table 1), corresponding to the horizontal dotted line. For ^^Ca, we have adopted 
the initial valne ^^Ca/^'^Ca = 3.9 x 10“^, more likely than the canonical ratio depicted in 
Table 1 (see text). 




Table 1. Initial Values of Radionuclides Considered in the Present Paper. 


R 

T. (Ma) 

D 

S 

R/S 

ref 

^Be 

53 d 

"Li 

^Be 

6.1 X 10-3 

a 

^^Ca 

0.1 

41k 

^oCa 

1.5 X 10-3 

b 

36C1 

0.3 

36g 

35C1 

1.6 X 10“^ 

c 

26 Al 

0.74 

26Mg 

2"A1 

4.5 X 10-5 

d 

i°Be 

1.5 

10b 

^Be 

1 X 10-3 

e 

eope 

1.5 

60Ni 

56Fe 

5-10 X 10-" 

f 

53Mn 

3.7 

53cr 

55 Mn 

4 X 10-5 

g 


Note. — R, D and S stand respectively for the ra¬ 
dionuclide of interest, its daughter and the stable refer¬ 
ence isotope. The R/S values are observed or inferred 
in CAIs. The original solar system values may some¬ 
times be higher (see text), (a) Chaussidon, Robert & 
McKeegan (2005). (b) Srinivasan et al. (1994). (c) 
Lin et al. (2005). (d) MacPherson, Zinner & Davis 
(1995). (e) McKeegan, Chaussidon & Robert (2000). 
(f) Tachibana et al. (2005); Mostefaoui, Lugmair & 
Hoppe (2005). (g) Birck & Allegre (1985). 







Table 2. Variation of the Beryllium-7 Irradiation Yields with the Spectral Parameters. 




Gore-mantle composition 

Ghondritic composition 



p = 2.7 

p = 4 

p = 5 

p = 2.7 

p = 4 

p = 5 

3He/iH = 

= 0 

1.0 (0) 

2.3(-l) 

9.6(-2) 

9.9(-l) 

2.2(-l) 

9.3(-2) 

^He/^H = 

= 0.01 

1.1 (0) 

2.6(-l) 

1.3(-1) 

1.0(0) 

2.5(-l) 

1.3(-1) 

^He/^H = 

= 0.1 

1.5(0) 

5.6(-l) 

4.5(-l) 

1.4(0) 

5.4(-l) 

4.3(-l) 

^He/iR = 

= 0.3 

2.4(0) 

1 . 2 ( 0 ) 

1.1(0) 

2.3(0) 

1.2(0) 

1.1(0) 

3He/^H = 

= 1 

5.7(0) 

3.5(0) 

3.6(0) 

5.5(0) 

3.4(0) 

3.5(0) 


Note. — The ^Be/®Be ratio is normalized to its measured meteoritic value, 
^Be/®Be = 6.1 x 10“^. Except for the spectral parameters, all other parame¬ 
ters (irradiation time, cosmic-ray flux, etc.) are kept the same as in Gounelle 
et al. (2001) (see section 2.1). The bold type number corresponds to the pre¬ 
ferred case of Gounelle et al. (2001) with p = 4 and ^He/^H = 0.3. 








